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( 1 ) 



Winninghoff [6], in his doctoral thesis investigated the effect on 
the process of geostrophic adjustment [2] in one dimension of several 
different spatial discretization schemes. His results were later in- 



corporated by Arakawa and Lamb [1] . The model he chose, which had 



been previously investigated by Rossby [4], Cahn [3] and others, was 



the linearized one-dimensional shallow water equations with no mean 



flow, in an infinite region: 




( 1 ) 




( 2 ) 




(3) 



where u is the perturbation velocity in the x direction, v is the 
perturbation velocity normal to the x direction, H and h are the mean 
and perturbed heights of the free surface, respectively, and g > 0 and 
f > 0 are gravitational and Coriolis parameters, respectively. The purpose 
of this note is to point out certain errors in Winninghoff T s solutions, 
present the correct solutions, and reexamine his conclusions in light of 
the correct solutions. 

Winninghoff followed earlier studies in solving (l)-(3) by eliminating 
between the equations to arrive at: 



then solving (4) by a Fourier Transform approach. After solving (4), 
solutions for h and v are obtained by substitution into (2) and (3) , 




(4) 



( 3 ) 



although closed-form solutions are not presented in some of the papers. 
Note the dispersive character of (4) is clearly seen by assuming a wave 
solution: 



u(x,t) = A e 



i (kx-vt) 



which implies: 



v 



= f z + gHk z = f z (1 + X 2 k 2 ), X = /iH/f . 



( 5 ) 



In [5], we presented a method where (1) — (3) were transformed directly, 
then the transformed problem was solved as a system of coupled ordinary 
differential equations?" If we denote transforms by an overhead tilde, e.g. 

CO 

u(k,t) = / u(x,t)e ''■^ x dx , 



then (1) — (3) reduce to 



= fv - ikg h 



d v 
dt 

dh 

dt 



= -fu 



= -ikH u 



( 6 ) 



Now, if we denote initial values by a subscript, e.g. 

co 

-ikx 



Uq ( k) = / u(x,0)e dx = u(k,0) , 



and use Winninghof f ' s initial conditions: 

u(x,0) = 0 
h(x,0) = 0 



v 



(V , -a < x < a 
(x,0) = l 0 

(o, otherwise , 



then using equation (12) of [5], the transformed solution is: 



( 7 ) 



( 4 ) 



u(k, t) = 



v(k, t) 
h(k,t) 



2fV n . i , \ 

0 . , . sin(ak) 

sin(vt) : 

v k 



ott ) k^gH , J sin(ak) 

2V„ < — j cos(vt)/ £ — - 

v v ) 



0 

2ikHfV 



Or-, / al sin(ak) 

2 — {1 - cos (vt ) } £ — - . 



v 

Winninghof f , using elimination, solved only for u(x,t) , 



( 8 ) 



uu 

u(x,t) = 2 ^ / u(k, t)e 



ikx 



dk 



(9) 



2faV A °° . , . . . , . 

0 j sin(ak) sinfvtl cos kx dk 

7T q ak v 



then solved for h(x,t) using the continuity equation. [If this is done 
to (9), or if the third equation in (8) is inverted directly, we find 
2HfV 0 00 n , u 

h(x,t) = / C - ^ S V sin(ak) sin(kx)dk . (10)] 

TT i. 2 

0 v 

The error in his solution occurs in converting from the differential case 
to the discretized case, using any one of the difference schemes he denotes 
as schemes A-E. (Figure 1) He performs this conversion by simply replacing 
v in (9) by its equivalent finite difference form, evaluating the integral 
numerically, then finding h(x,t) from the continuity equation. We now 
show this approach neglects the effects of spatial averaging of u(x,t) 
and v(x,t) that occurs in Schemes C and D. For example, in Scheme C, 
the finite difference form of (l)-(3) is: 



3u _ f [v(x+d/2,t) + v(x-d/2,t)] 
3t 2 



[h(x+d/2,t) - h(x-d/2,t)] 
d 



3v n [u (x+d/ 2 , t) + u(x-d/2,t)] 
3t 2 



(ID 



_ tt [u(x+d/ 2 , t) - u (x-d/2 , t) ] 
_ _ _ H _ 



( 5 ) 



If we transform this system, distinguishing from the differential 
system by using a subscript c on the variables, we have: 

dU c r . sin (kd/2) r 

— = f cos (kd/2) v c - ig (572) h c 

d\/ ^ 

= -f cos(kd/2)u c 

* h c sin(kd/2) ~ 

dE lH (d/2) "c ' 



( 12 ) 



This system is identical to (6), except that f is replaced by 
f cos(kd/2) and k is replaced by sin (kd/2) / (d/2) . Thus the 
solutions for u^(x,t) and h c (x,t) , corresponding to Arakawa’s 
initial condition, can be obtained by making a similar replacement in 
(12) of [5], to yield: 



2faV 



u Q (x,t) = 



0 



7T q ak 
00 (1-cos (v t)} 



sin(ak) s ^- n 



2Hf aV 



h c (t) = 



0 



cos(kx) cos(kd/2)dk , 



(13) 



/ 

0 



C 

sin(ak) 

ak 



sin (kx) s i n (kd/2) cos (kd/2) dk 



where 



V 2 = f 2 (1 + 4 (A/d) 2 sin 2 (kd/2) ) . 
c 



(15) 



Observe that (13) is not obtainable from (9) simply by replacing v by 

v . This is because of the cos(kd/2) term that arises from the 
c _ 

spatial averaging of v and u that occurs in the first and second 
equations of (11), respectively. Repeating these calculations for the 
other schemes will show that u(x,t) in Scheme D also requires a 
cos(kd/2) term, and that h(x,t) is given, in general, by 



2Hf aV 



h(x , t) = 



0 



/ D(k) 



(1 - cos vt) sin(ak) 
2 ak 



sin(kx)dk 



(16) 



where 



(14) 



( 6 ) 



Differential 



D(k) = 



and 



v>/ f = 



sin (kd) 
d 

sin (kd/2) 

(d/2) 

cos (kd/2) sin (kd/2) 
(d/2) 

cos ^ (kd/2) sin (kd/2) 
(d/2) 



sin(kd/ /2) 
(d//2) 



2 2 

[1 + AV] 



1/2 



[1 4 (A/d)^ sin^ (kd) ] 



1/2 



[1 4- 4(A/d)^ sin^ (kd/2)] 



1/2 



, Scheme A 
, Scheme B 
, Scheme C 
, Scheme D 
, Scheme E 



[cos^ (kd/2) 4- 4(A/d)^ sin^ (kd/2)] 



1/2 



[cos^ (kd/2) 4* (A/d)^ sin^ (kd) ] 



1/2 



[1 4- 2(A/d)^ sin^ (kd//2)] 



1/2 



, Differential 



Scheme A 



Scheme B 



Scheme C 



Scheme D 



, Scheme E 



In Figure 2, the integrals given by (16) are used to compute (h/H) at 
t = 80 hours in the same manner as done by Winninghoff, i.e., using 
Simpson’s rule with 600 subintervals in the interval (0,7r/d) . The 
relevant parameter values used were 

f = 10 ^ sec 1 



g = 10 m sec 



-2 



H = 10 m 



a = d = A/2 . 



( 17 ) 



(18) 



( 7 ) 



These graphs too should be compared to Figure 3, which shows the 
equivalent evaluation of (10), with only v(k) replaced by its equiv- 
alent finite difference value. (Figure 3 appears identical to Figure 
3.6 in Winninghof f . ) Note especially the effect on the solution to 
Scheme D. We conclude that Winninghoff obtained such ill-behaved 
solutions in this case because he neglected the cos(kd/2) term in 
the numerator of u(x,t) in these cases, and this term helps to 
control the singular behavior of [v(k)] ^ as k-vrr in Scheme D. 

Winninghoff eliminated Scheme D from further study due to the 
extremely noisy behavior shown by Figure 3. Since his calculation was 
erroneous, it can not be discarded for this reason, however, as he 
showed, it is still seriously deficient in terms of phase speed errors, 
and should probably be discarded for this reason. 
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FOOTNOTE 



1. A sgn term was erroneously omitted from certain expressions in 
[5], Thus, in [5], equation (13) should read, 

hg(x) = h(x, 0 ) ~~ / sgn(x - s)e S |^(s, 0 ) - v(s, 0 )}ds , 

2\ f -oo 



the equation immediately above (28) should read 

“$|x| A 



i- / 

2 ir J 



ic 



2 2 

1 + Ac 



ikx ,, - Ax e 

e dk = - 



sgn(x) 



A [A 2 + (Ax ) 2 ] 1/2 n=-°° 



l «S(x - ( 2n - 1) Ax) , 



equation (30) should read 



h (x, 0 ) = h(x , 0 ) + ~ 5 — yTo 

b 2Af [A Z + (Ax Z )] ±/Z 

x \ {sgn((2n - l)Ax)e ^ ‘ l)Ax| X ^ _ ( 2 n - 1) Ax ,0) (2Ax) } , 

n— —oo 



and equation ( 2 . 8 ) should be identical to the second equation in this 
note. 
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